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Abstract 

We discuss the question of inverse symmetry breaking at non-zero temperature 
using the exact renormahzation group. We study a two-scalar theory and concentrate 
on the nature of the phase transition during which the symmetry is broken. We also 
examine the persistence of symmetry breaking at temperatures higher than the critical 
one. 
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It has been known for a long time 0, @, 0], that simple multiscalar models can exhibit 
an anti-intuitive behaviour associated with more broken symmetry as the temperature is 
increased. We refer to this behaviour as inverse symmetry breaking. This possibility may 
have remarkable consequences for cosmology, solving the problem of topological defects, 
thanks to the fact that the phase transition leading to their formation may never have 
occurred during the thermal history of the Universe [§. The existence of the phenomenon 
for supersymmetric theories has been discussed in refs. 0. However, doubts have been 
raised on the validity of these results, which are based on the one-loop approximation to the 
non- zero-temperature effective potential. This approximation is known to be unreliable for 
the discussion of many aspects of phase transitions. Recently the effect of next-to-leading- 
order contributions within perturbation theory has been investigated in ref. [^]. This has 
been done through the study of gap equations, which are equivalent to a resummation of 
the super-daisy diagrams of the perturbative series. Large subleading corrections have been 
identified, which lead to a sizeable reduction of the parameter space where inverse symmetry 
breaking occurs. The question of inverse symmetry breaking has also been studied through 
the use of the renormalization group, with similar conclusions 0. A variational approach 
has been employed in ref. [§. Contrary to the results of the above studies, a large- 

analysis seems to indicate that symmetry is always restored at high temperature ||^. 
However, the validity of this claim has recently been questioned in ref. A finite- lattice 



calculation also supports the symmetry restoration at sufficiently high temperature [p 
even though the relevance of this result for the continuum limit is not clear. 

In this letter we study the question of inverse symmetry breaking by employing the 
exact renormalization group [l^. Our formalism is similar to that of ref. [0. Our ap- 



proach, however, is based on the real-time formulation of high-temperature field theories. 
We investigate a two-scalar theory and we identify the universal behaviour associated with 
the symmetry-breaking phase transition. Our study is based on an evolution equation for 
the potential of the non-zero temperature theory. Using a polynomial ansatz for the po- 
tential, we solve this equation and verify the conclusion of ref. that inverse symmetry 
breaking is confirmed by the renormalizat ion-group approach. Moreover, we explore the 
parameter space that leads to inverse symmetry breaking and compare it with the pertur- 
bative predictions. In order to study the phase transition, we go beyond the polynomial 
ansatz and consider a general dependence of the potential on the field that develops an 
expectation value when the symmetry is broken. The resulting partial differential equation 
is solved numerically with the use of appropriate algorithms |]13|. The Wilson-Fisher fixed 
point of the effective three-dimensional theory is shown to govern the dynamics near the 
second-order phase transition. 

We consider the simplest model that exhibits inverse symmetry breaking: a two-scalar 
model with Z2 x Z2 symmetry. The tree-level potential is given by 

VtMi. 02) = + + \^x4>\ + ^A20^ - ^Ai20?0^ . (1) 

This potential is bounded for Ai,2 > and 

A1A2 > A?2 . (2) 
The thermal correction to the above potential at the one-loop level is given by the well- 
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known expression 
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Here J^f^ indicates the matrix of second derivatives of the tree-level effective potential 
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When both eigenvalues of Ai^j. are much smaller than (which happens for sufficiently 
small couplings), the leading field-dependent correction takes the form 

rp2 
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3Ai - A12 < 
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which can be consistent with the stability condition of eq. (|]), the thermal correction for 
the mass term of the (pi field is negative. If the system is in the symmetric phase at zero 
temperature with > 0; there will be a critical temperature T^^ = 12m^/(Ai2 — 3Ai) 
above which the symmetry will be broken If the system is in the broken phase at T = 0, 
the symmetry will never be restored by thermal corrections. 

Our aim is to discuss the above scenario in the context of the Wilson approach to the 
renormalization group. The main ingredient in this approach is an exact flow equation 
that describes how the effective action of the system evolves as the ultraviolet cutoff is 
lowered. We consider the lowest order in a derivative expansion of the effective action, 
which contains a general effective potential and a standard kinetic term. At non-zero 
temperature this approach can be formulated either in the imaginary-time [0 or in the 
real-time formalism |jl6|. In the latter formulation, the evolution of the potential lowering 



the cutoff scale A is given by the partial differential equation |T6| 
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Notice the formal similarity with the Dolan-Jackiw one- loop result of eq. (|^). The main 
difference is that the 'running' mass matrix in the exponential of eq. replaces the tree- 
level one in eq. (0). The initial condition for the above equation, at a scale Aq ^ T, is the 
renormalized effective potential at zero temperature. In this work we consider small quartic 
couplings, so that the logarithmic corrections of the zero-temperature theory can be safely 
neglected. The initial condition for the evolution is a zero-temperature potential given 
by eq. (|1]). Integrating the evolution equation (^, we obtain the non-zero-temperature 

^ For 3A2 — A12 < the symmetry is broken in the (f>2 direction. The discussion of this case is completely 
analogous to the one we consider. Notice that it is not possible to break the symmetry in both the 0i and 
02 directions, because of the stability condition of eq. (0). 
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effective potential in the limit A — >^ 0. In the approximation that the 'running' mass matrix 
on the r.h.s. of the evolution equation is taken to be constant, equal to the tree- level 
mass matrix, the integration reproduces the perturbative result of eq. The non-trivial 
behaviour that we describe in the following paragraphs is obtained when the full scale- 
dependence of the mass matrix is taken into account. An interesting limit corresponds to 
the evolution at scales A <^ T. Let us denote by Ai\ ^ the eigenvalues of the mass matrix 

that satisfy Ai\ ^ ^ A^. The remaining eigenvalues correspond to decoupled massive 
modes that do not contribute to the evolution at the scale A fl^. Keeping the leading 



contribution in the r.h.s. of eq. and omitting the field- independent terms, we obtain 
d A3 
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The rescaled potential Va/T has dimensions (mass)^ and its evolution is typical of that 
of a three-dimensional theory |18|. In the limit A -C T, dimensional reduction takes 
place, and a four-dimensional theory at non-zero temperature behaves as an effective three- 
dimensional one at T = 0. We can cast eq. (||) in a form that does not depend explicitly 
on the scale A by defining the dimensionless parameters: 
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The dimensionless mass matrix Ml related to the rescaled potential V through an 
equation analogous to eq. (H). The evolution equation now reads 



3V^-^Tr{log(l+[A<i]J} , 



where again we have omitted the field-independent terms. The scale-invariant (fixed-point) 
solutions of the effective three-dimensional theory can be obtained from the above equation 
for KdVK/dk = 0. 



Finding the solution of eq. (^ is a difficult task. An approximate solution can be 
obtained ||T5| , 0] by expanding the potential in a power series in the fields. In this way the 
partial differential equation (|^) is transformed into an infinite system of ordinary differential 
equations for the coefficients of the expansion. This system can be solved approximately 
by truncation at a finite number of equations. In effect, the potential is approximated by 
a finite-order polynomial. As a first step, we follow this procedure and define the running 
masses and couplings at the origin 
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The corresponding evolution equations can be obtained by differentiating eq. 
neglecting the higher derivatives of the potential. We find 
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with 



A3 Ar(cUi,2) ^ 9Ci2 2 a2 2 .X 

C,,2 = —.,^-^, D,,2 = Tr-^, u;'=A' + ml,, (14) 
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and N{uj) = [exp(co'/T) — 1] ""^ the Bose-Einstein distribution function. For uji^2 <^ T we 
have 

A3 T 

Cl,2 T^T9~, 2~ ' (-'-^) 

47r^ y\2 _l_ ^^^^ 

and the above equations agree with those considered in ref. in the same hmit. For 
^1,2 ^ T there is no running, because of the exponential suppression in the Bose-Einstein 
function. 



We have solved numerically the system of equations ([I3D and determined the range of 
zero-temperature parameters that lead to inverse symmetry breaking. In fig. 1 we present 
the results for a zero-temperature theory with positive mass terms mf(Ao) = mKAo) 
and A2(Ao) = 0.3. The temperature has been chosen much higher than the critical one 
(T = 500mi(Ao)). The system (|13|) has been integrated from Aq ^ T down to A = 0, 
where the thermally corrected masses and couplings at non-zero temperature have been 
obtained. A negative value for the mass term ml at A = has been considered as the 
signal of inverse symmetry breaking. This has been achieved for the parameter range of 
Ai(Ao), Ai2(Ao) above the line (a) in fig. 1. In the same figure we plot the stability bound 
of eq. (H). The allowed range is below the line (b). We also include the perturbative 
prediction of eq. @ for the range that leads to inverse symmetry breaking. It lies above 
the line (c). The phenomenon of inverse symmetry breaking is confirmed by our study, in 
agreement with ref. [0, where the imaginary-time formulation of the renormalization-group 
approach has been used. We observe that the renormalization-group treatment eliminates 
a large part of the parameter space allowed by perturbative theory. This is in agreement 
with the results of ref. 0, where the gap-equation approach has been followed. Notice that 
lines (a) and (c) approach each other near the origin, where perturbation theory becomes 
more reliable. 

The reliability of our conclusions crucially depends on whether the solution of the system 
of truncated equations ([T3| ) provides an approximate solution to the full partial differential 
equation (^. In ref. 0, this has been checked by increasing the level of truncations and 
verifying the convergence of the results. We follow here a different approach that relies on 
the numerical integration of eq. through the algorithms discussed in ref. These 
algorithms have been used for the integration of the evolution equations for potentials that 



depend on one field only [13, IS]. The generalization to the two-field case is straightforward. 
However, limitations in computer time have prevented us from reaching full numerical 
stability for the results. For this reason, we restrict our discussion of eq. along the 0i 
axis, which is the direction of expected symmetry breaking for our choice of couplings. We 
approximate the potential by the expression 

^a(0i, = ^a(0i) + lrnl{A)<Pl + ^\2{A)<Pi - ^Ai2(A)0?02 _ (^g) 

The evolution of ml{A), A2(A) and Ai2(A) is determined through the truncated eqs. (JTB)). 
However, the full 0i dependence is preserved through the numerical integration of eq. (|^, 



4 



with the eigenvalues of the mass matrix Ai\ given by 



9Va(0i) 



Mi\^ = ^^ and [Mi\^ = mi{A)-X,,{A)<Pi. (17) 

This treatment permits a rehable study of the order of the symmetry-breaking phase tran- 
sition. The appearance of secondary minima of the potential at some point in the evolution 
can be studied in detail. As a result, we can distinguish between first-and second-order 
phase transitions. This is not possible when local expansions of the potential, such as the 
one leading to eqs. (|TB|), are employed. 



In fig. 2 we present the evolution of the potential for zero-temperature parameters 
ml{Ao) = m2(Ao), Ai(Ao) = 0.01, A2(Ao) = 0.3, Ai2(Ao) = 0.05. The location of this 
parameter choice on the plot of fig. 1 is denoted by a black square. It is within the region 
for which inverse symmetry breaking is expected. The temperature is chosen very close to 
the critical one: Tcr/mi{Ao) ~ 33.3. This value is in very good agreement with the result 
of ref. 0], for the same choice of zero-temperature parameters. It deviates significantly 
from the perturbative prediction Tcr/rriilAQ) ~ 24.5. We use the rescaled variables defined 
in eqs. (p!oD, which permit the identification of the fixed points that may be relevant for 
the phase transition. We plot the derivative dV\/dpi of the potential as a function of the 
variable pi = 0i/2, for decreasing A. The early stages of the evolution of the potential, 
when its curvature at the origin becomes negative, are not clearly visible in this plot. The 
reason is that dVty/dpi is very small for large A. The important point in this figure is the 
flow of the potential towards a scale-invariant solution (marked by WF) during the later 
stages of the evolution. This solution corresponds to the Wilson-Fisher fixed point of the 
effective three-dimensional theory. For A — > the theory leaves the fixed point and flows 
towards the phase with symmetry breaking. The evolution can also be described in terms 
of the parameters defined in eqs. (|T2p. Their rescaled versions 



"^l2(A) = '"^^2"' ^l,2,12(A) ^ y^^j 

are obtained from the rescaled potential of eqs. (|^) through relations analogous to eqs. (|T^ . 
(Notice that fn{{A) is given by the value of dV\/dpi for pi = in fig. 2.) We present the 
evolution of all these parameters in figs. 3 and 4. In both these plots we observe that 
the mass term m\{A) turns negative at some scale A ~ T. The mass term m\{A) stays 
positive and grows as A decreases. For A <C T we observe that m\{A) ^ A^ and the 02 
field decouples. Eventually, for A ^ 0, the renormalized parameters m\{A) and A2(A) take 
constant values. After the decoupling of the 02 field, only the fluctuations of the (pi field 
contribute to the evolution. In this regime, the evolution equation for the potential can be 
written in the form of eq. (|Tll), with only the eigenvalue }A\ ^ = c}^Va/90^ contributing. 
The scale-invariant solution of the resulting equation corresponds to the Wilson-Fisher fixed 
point of the one-scalar three-dimensional theory. Its explicit form can be seen in fig. 2. We 
conclude that the role of the 02 field is to trigger the inverse symmetry breaking by inducing 
a negative mass term (A) at the early stages of the evolution. It subsequently decouples 
and the dynamics of the phase transition is governed by the Wilson-Fisher fixed point of 
the one-scalar three-dimensional theory. The resulting phase transition is of second order. 
Its universal behaviour can be parametrized in terms of critical exponents and amplitudes. 



(^) X ..^_Aw2(A)T 



For a detailed discussion see ref. [15 
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The evolution close to the fixed point is apparent in fig. 3, in the range where the rescaled 
parameters m^(A), Ai(A) become constant. Their fixed-point values are ml jp ~ —0.44, 
Xijp — 3.2. In this regime, the parameters m^(A), Ai(A) evolve towards zero according 
to m^(A) = rhfjpA'^, Ai(A) = AijpA/T (see fig. 4). This explains why the fixed point 
is relevant very close to the phase transition. If the temperature is chosen such that the 
evolution stays close to the fixed point for a long 'time' t = log(A/T), the curvature of the 
potential at the origin is very small at the end of the evolution. Through sufficient fine- 
tuning of the temperature, the curvature can be made arbitrarily small. This is the criterion 
for the occurrence of a second-order phase transition. Notice that the renormalized quartic 
coupling Ai(A = 0) is also expected to be zero at the critical temperature fl^. According 
to eqs. ( [T3| ) the coupling Ai2(A) evolves as A12 ~ A", with a = 3Aijp/27r^(l + m^jp)^ ^ 1.6. 

For A2(A) we obtain A2(A) = A2(Ai) + c (a^^-i - A-''"^) /(2a - 1), where A^ is the scale 
at which the fixed-point solution is initially approached, and c a calculable constant. The 
predicted evolution of Ai2(A) and A2(A) is confirmed by fig. 4. If we had used the truncated 
equations (|I^ for the discussion of the evolution of m^(A) and Ai(A), as in ref. 0, we 
would have obtained a = 1/3. This would have predicted an unphysical, singular behaviour 
A2(A) ~ — A~^/^ for A — + 0. Our solution of the evolution equation for the full potential 
along the (pi axis resolves this problem, and describes the decoupling of the 02 field properly. 

In conclusion, we have used the real-time formulation of the exact renormalization group 
in order to study the question of inverse symmetry breaking in the context of the two-scalar 
theory of eq. (|I|) . We have verified the conclusion of ref. [0] that inverse symmetry breaking 
is confirmed by the renormalizat ion-group approach. We have also determined the param- 
eter range that leads to inverse symmetry breaking, through appropriate truncations of the 
partial differential equation that describes the evolution of the potential. This parameter 
range is significantly smaller than the perturbative prediction, in agreement with ref. 
We have also obtained a numerical solution of the evolution equation for the full poten- 
tial along the 0i axis, without relying on truncations in that direction. In this way we 
have obtained a detailed picture of the symmetry-breaking phase transition. The 02 field 
triggers the symmetry breaking along the 0i direction and subsequently decouples. The 
phase transition is of second-order, governed by the Wilson-Fisher fixed point of the effec- 
tive three-dimensional theory. Our improved treatment gives no indication of symmetry 
restoration for the range of temperatures above the critical one, in which we obtain stable 
numerical solutions. 
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Figures 

Fig. 1 The parameter space that leads to inverse symmetry breaking (A2 = 0.3). 

Fig. 2 The derivative = dV\/ dpi of the potential as a function of pi = (pl/2 for decreasing 
A, for a theory very close to the phase transition. 

Fig. 3 The evolution of the rescaled masses and couplings for a theory very close to the 
phase transition. 

Fig. 4 The evolution of the masses and couplings for a theory very close to the phase tran- 
sition. 
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